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1. Introduction 

In the past few years our understanding of QCD at finite temperature has advanced consider- 
ably, due to improvements in the lattice discretizations (use of improved actions), improvements 
in the simulation algorithms, and increases in computational resources. A sketch of our present 
understanding of the QCD phase diagram in the m uc i — m s plane is shown in Fig. |l|. 

N f = 2 Pure 




Figure 1: Sketch of the QCD phase diagram in the m m \ — m s plane. 

In this talk I review recent progress in the determination of the QCD phase diagram at finite 
temperature and zero baryon density, in investigations of the nature of the transition or crossover 
from the hadronic phase at low temperatures to the quark-gluon plasma (QGP) phase at high tem- 
peratures and in the determination of the equation of state (EOS). For earlier recent reviews see 
Refs. [jl]], while a summary of other properties of the high-temperature phase is given in the talk by 
T. Hatsuda [||. 

2. Simulation choices and improvements 

As in any numerical simulation, for finite temperature studies one must makes choices of 
the action and parameters to be simulated. The most important choice, for dynamical fermion 
simulations, is the choice of fermion action, each of which has advantages and drawbacks. To 
weigh those we recall that the finite temperature transition or crossover is driven, for small quark 
masses, by the restoration of chiral symmetry in the high temperature phase. 

(i) Wilson-type fermions, including clover fermions, have the advantage that they are local (in fact 
"ultralocal") for any number of fermions, while the fermion determinant is positive for even num- 
bers of fermions. Their big disadvantage is that the chiral symmetry is explicitly broken and the 
chiral limit therefore not protected. This makes the study of chiral symmetry restoration cumber- 
some and difficult. It becomes really meaningful only in the continuum limit. 

(ii) Staggered fermions have the advantage that they have a (remnant) chiral symmetry. The chiral 
limit is thus protected. In addition, they are comparatively cheap to simulate. The main disadvan- 
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tage of staggered fermions is the need to use the "fourth root trick" when the desired number of 
fermions of a given mass is not a multiple of four. 

(iii) Overlap and domain wall fermions (at least for sufficiently large fifth dimension L s ) have the 
advantage of good chiral symmetry and the protection of the chiral limit, and of being local for 
any number of fermions (at least for sufficiently small lattice spacing a). Furthermore, for overlap 
fermions, the fermion determinant is positive for any number of fermions [§]. The disadvantage of 
these chiral fermion discretizations is that they are expensive to simulate. This problem is acerbated 
by the fact that one might need lattices with N t > 8 to have lattice spacings small enough in the 
transition/crossover region for the fermion action to be sufficiently local. In judging this drawback, 
one should keep in mind that, for the computation of the EOS, for example, the scaling of the 
costs of a simulation is worse, by a factor ~ aT A , than for typical zero temperature simulations, 
because the observable, obtained from the difference of a finite and a zero temperature simulation, 
decreases, in lattice units, as a 4 , while the error decreases much slower. 
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Figure 3: The pion mass splitting, (m^ — 
m 2 n )/T 2 as function of rn^/T 2 at the lattice spac- 
ing corresponding to the finite temperature transi- 
tion/crossover point, from Ref. [H], 



Most of the recent dynamical fermion thermodynamics simulations have used improved ver- 
sions of staggered fermions. The improvement aims to reduce taste symmetry breaking, by using 
some version of "smeared" or "fattened" gauge links in the nearest neighbor hopping term. This 
is done for all three versions currently under investigation, p4 (Bielefeld, RBC-Bielefeld) [Q, asq- 
tad (MILC) ^ and stout-link (Wuppertal-Budapest) [^]. The first two also aim to improve the 
dispersion relation, which implies improved behavior of thermodynamic quantities in the high 
temperature limit {i.e., for free fermions) by including three-hop terms: straight, the Naik term, for 
asqtad, and "knight moves" for p4. For free fermions, and hence in the high temperature limit, the 
link fattening becomes inoperative. Thus stout-link fermions act like standard staggered fermions, 
and asqtad fermions like Naik fermions. The two forms of improvement are illustrated in Figs. || 
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and|. 

The other major advance, in the last couple of years, has come from the algorithmic side. 
Because of the fractional power of the fermion determinant, needed for the fourth-root trick, the 
usual, exact HMC algorithm [j|] can not be used. Instead, the inexact hybrid R algorithm [Q], with 
stepsize errors of 0{e 2 ), was employed. This obstacle was recently overcome with the invention 
of the exact RHMC algorithm [0. As discussed in the talk by M. Clark [jlT]] the RHMC algorithm 
is not only exact but allows various other improvements - multiple time-step integration schemes, 
multiple pseudofermion fields, etc. - that can speed up the simulations by factors of 2 up to 8 
compared to the R algorithm. 

A couple of exploratory studies with Wilson-type fermions were presented at this conference, 
using a hypercube action [12] and with twisted-mass Wilson fermions [13]. But in the rest of this 
talk I will concentrate on simulations with (improved) staggered fermions. For the purpose of this 
talk the validity of the fourth-root trick, discussed in detail by S. Sharpe [ 14], will be assumed. 



3. The phase diagram 

3.1 The physical point 

Determining the nature of the finite temperature transition or crossover at the physical point is 
not only of interest in its own right, but has implications on the possible phase diagram with finite 
chemical potential (see e.g. the talks by C. Schmidt [ |l5| ] and M. Stephanov[[l^]). The Wuppertal- 
Budapest group [|l7|] made a systematic investigation using stout-link improved fermions and the 
exact RHMC algorithm. They used a physical strange quark mass and physical (degenerate) light 
quark masses. The chiral susceptibility for N t = 4 and 6 and various volumes is shown in Fig. |j. 
It shows no sign of increasing with volume as would be expected for a genuine phase transition, 
indicating existence of only a crossover. 
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Figure 4: The chiral susceptibility for Nt =4 (left) and Nt =6 (right) and various volumes, from Ref. [|l7|]. 

A definite statement, however, needs an extrapolation to the continuum and infinite volume 
limits. For a meaningful extrapolation renormalized observables have to be considered. The 
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Wuppertal-Budapest group does this by subtracting the T = value, Ax = %{T) — %{T = 0), 
to cancel a potential additive divergence, and then multiplies with m 2 q to obtain an RG invariant 
observable. This is then extrapolated to the continuum limit for fixed aspect ratio N s /N t = 3,4 and 
5, see Fig. ||. 




Figure 5: The continuum extrapolation of T 4 / (rn^Ax) for fixed aspect ratio N s /Nt = 3,4 and 5, from 
Ref. [Q. 

The resulting continuum renormalized susceptibility stays finite and non-zero in the infinite 
volume limit. It would diverge for a genuine phase transition. Hence we now have convincing 
evidence that, in nature, QCD has a crossover at finite temperature and zero baryon chemical po- 
tential. 

3.2 The second order boundary line 

De Forcrand and Philipsen mapped out the second order critical line that separates the first 
order region at small quark masses from the crossover region at intermediate quark masses in the 



T^ud — m s plane of the phase diagram [18]. They used standard, i.e., unimproved, staggered fermions 
and the exact RHMC algorithm. 

At fixed am s they determined the location of the second order point, as function of am uc i, by 
requiring that the Binder cummulant, B^(m'^) attain its critical Ising value, ~ 1.6. The locations, 
with errors, are shown in Fig. |[ The arrows in the plot indicate parameters values at which zero- 
temperature simulations were performed to determine meson masses and lattice spacing. De For- 
crand and Philipsen found 

• For the physical strange quark mass, the second order boundary occurs at non-vanishing 
m< ud < m ud S > l '- e > tne P n y sica l point is in the crossover region, in agreement with the result 
of the Wuppertal-Budapest group. 

2/5 

• The critical line is consistent with m c /"(m uc i) = m'" — cm ud , i.e., the behavior that is ex- 
pected if a tricritical point exists. 

One should emphasize that these results are most likely qualitative, only. Unimproved stag- 
gered fermions were used on lattices with N t = 4, were the lattice spacing in the crossover/transition 
region is quite large, a ~ 0.27 fm. Therefore lattice effects could be significant. For exam- 
ple, at the second order critical point for three degenerate flavors, de Forcrand and Philipsen 
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find m n /m p = 0.304(2), whereas the Bielefeld group, using p4 fermions, concluded that [19] 



m K /m p < 0.18, the physical value of the mass ratio. For additional results with three degener- 



ate flavors see also the contribution by Cheng Q20|]. 

3.3 Massless 2 flavor #QCD 

Kogut and Sinclair continued their study of massless Nf = 2 ^QCD on N t = 8 lattices with 



N s = 12, 16 and 24 [glj. #QCD contains an irrelevant chiral 4-fermion interaction, which makes 
the fermion matrix non-singular in the limit of zero (bare) quark mass while preserving the chiral 
symmetry of the Lagrangian, thus enabling the study of spontaneous chiral symmetry breaking in 
the massless limit. For numerical simulations, the 4-fermion interaction is made quadratic in the 
fermions by introducing an auxiliary (o,n) field. 

Kogut and Sinclair used (standard) staggered fermions which, at finite lattice spacing have a 
U(l) x U(l) chiral symmetry, expected to be spontaneously broken to U(l). Hence, if the chiral 
symmetry restoration transition is second order, one expects the transition to be in the 3-d 0(2) 
universality class. Therefore, one wants to compare to the 3-d 0(2) spin model. However, the 
magnetization of the 0(2) spin model has large finite size effects, as seen in Fig. (left). Hence, 
instead of comparing to infinite volume 0(2) behavior, Kogut and Sinclair compare to finite volume 
behavior for a "best matched size". The strategy, thus, is to find the volume for which the 0(2) 
magnetization gives the best fit 

{\j/y(p))=b{M(a(P-p c ) + T c )) i.e. 1/7 = T = a{fi + T C . (3.1) 

of the chiral condensate, for a given QCD lattice size. 

Such fits work well, see Fig. ^ (right), strongly suggesting that the #QCD data are compatible 
with the 0(2) spin model. The corresponding 0(2) volumes are small, 6 3 for 16 3 x 8 and 8 3 for 
24 3 x 8. Thus, one should not attempt comparisons with large volume critical behavior. Kogut 
and Sinclair find, in retrospect, that such fits also work for their earlier N t = 6 data. Hence the 



conclusions of their earlier work, Ref. [22], should be modified accordingly. 
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0(2) spin model L 3 LATTICE 
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Figure 7: The magnetization of the 0(2) spin model in the transition region for various volumes (left) and 
an example of a best volume fit of (yy) to the 0(2) spin model magnetization (right), from Ref. [EJJ. 



3.4 2 flavor LQCD with KS quarks at N, = 4 



The Pisa group [23] continued their investigation of the phase transition of Nf = 2 QCD with 
standard staggered quarks at N t = 4. To check their previous results \&4], obtained with the inexact 
R-algorithm they performed comparisons with the exact RHMC algorithm. They found that the 
systematic ff(e 2 ) step-size errors are comparable to the statistical ones. 





r Order 



Figure 8: Test of first order finite volume scaling with miL kept constant, from Ref. []23[]. 



They then tested their hypothesis of a first order transition with a finite size scaling study of 
the specific heat, Cy, and chiral susceptibility, % m , with the scaling variable miD' h kept constant for 
the choice yj, = 3, the value for a first order transition, see Fig. || The specific heat scales nicely, 
but not the chiral susceptibility. They speculate that the reason for the latter might be the large mass 
needed for the smaller volume, am^ = \^ = 0.1068. 

On the other hand, T. Mendes [^] recently compared previous data by the MILC collabora- 
tion [[26]] and her own new data at a heavier mass, am = 0.075, with the (infinite volume) 0(4) 
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Figure 9: Comparison of Nf = 2 QCD data with staggered fermions atN t =4 to the 0(4) scaling function, 
fromRef. & 



scaling function, as shown in Fig. ||. Mendes inferred that the scaling works quite well, especially 
at the heavier mass. 

I would conclude that the issue of the order of the phase transition for 2-flavor QCD, even with 
staggered fermions at N t = 4, and certainly in the continuum limit, remains an open question. 

3.5 T c for full QCD 



The RBC-Bielefeld collaboration [ |27| , [28| ] has recently performed a systematic study of the 
crossover temperature in full QCD, i.e., physical strange quark mass, and degenerate light up and 
down quark masses, m/. They used p4 fermions and the exact RHMC algorithm for simulations 
with N t =4 and 6, several light quark masses mi, and fairly large volumes, 2 < N s /N t < 4. 

For each simulation parameter set, they located the crossover point from the peak position 
of the chiral and Polyakov loop susceptibilities. The two peak positions appeared to coincide at 
large volume, suggesting that the crossovers for deconfinement and chiral symmetry restoration 
occur simultaneously. So the authors averaged over the two determinations at finite volume. T c , in 



units of r Q 1 is shown in Fig. 10 as function of mpsrQ. The data at finite lattice spacing, a« 1 /N t , 



and unphysical light quark masses are then extrapolated to the continuum limit and physical point, 
given by nipsro = 0.321(5) — ro was taken from Ref. Q29Q: ro = 0.469(7) fm — using the form 



mi ,N t 



T c r +A(m PS r ) d + B/N? 



(3.2) 



Here, d = 1.08 would be expected for a second order phase transition (in the 0(4) universality 
class) and d = 2, i.e., a linear dependence on mi, for a first order transition. At the physical point 
they find T c r$ = 0.457(7)^3, with the second error coming from using d = 2 and d = 1 in the 
extrapolation. 

In physical units, the transition/crossover temperature is T c = 192(7) (4) MeV . This is about 
12% larger than the value obtained by the MILC collaboration, using asqtad fermions, the R- 
algorithm, lattices with ,/V, = 4, 6, 8 and N s /N t = 2, and a combined chiral/continuum limit extrap- 
olation as in eq. (Q. The MILC collaboration obtained T c = 169(10) (4) MeV []7p. While they 



worked in units of the scale r\ , the conversion to physical units was done with compatible values 
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Figure 10: T c in units of r Q 1 as function of mpsro, 



from Ref. [27 1. The vertical line shows the location 
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Figure 11: Recent lattice determinations of T c , 



of the physical value. 



collected by R Petreczky in Ref. [ 30 1 



and does not contribute to the disagreement. A collection of recent determinations of T c , assembled 



by R Petreczky in Ref. [ 30], is shown in Fig. 1 1 



4. The equation of state 

After establishing the phase diagram, including the order and nature of the phase transitions 
and/or crossovers, one would like to understand the nature of the different phases. One of the 
basic quantities for this is the equation of state (EOS), the pressure p(T), the entropy density 
s(T) = dp(T)/dT and the energy density s(T) = Ts(T) — p(T). Besides its intrinsic interest as 
a fundamental property of QCD, the EOS is of phenomenological interest. For example, it is an 
import input in hydrodynamical models of the QGP, often used to try to interpret results from 
heavy-ion collision experiments, such as the observed elliptic flow. For a quantitative understand- 
ing, obviously, a quantitative understanding of the EOS is necessary. 

4.1 Low temperature behavior 
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Figure 12: Comparison of the EOS from lattice simulations with a hadron resonance gas model with masses 



adjusted to the lattice quark masses used, from Ref. [31 1. 
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In the low temperature hadronic phase, the excitations are weakly interacting hadrons, includ- 
ing all resonances. This suggests that a hadron resonance gas model (HRG) should give a good de- 
scription of the EOS. This was recently tested by a comparison with lattice data using p4-fermions 
and N t = 4 Since the light quark mass used was larger than the physical light quark mass, the 
hadron masses of the HRG where adjusted to match those on the lattice. With this adjustment the 



HRG works surprisingly well, even around the crossover temperature, as can be seen in Fig. [12. 

4.2 The EOS at high temperatures 

At (very) high temperatures the running coupling g(T) vanishes (logarithmically in T), and 
one expects that the EOS should be computable in perturbation theory. PT, though, is afflicted with 
infrared divergences, requiring resummation of certain classes of diagrams, giving e.g. an &(g 3 ) 
and an 0(g 5 ) contribution. At &{g^) the IR divergences cause all orders of PT to contribute, so 
that this order needs to be computed non-perturbatively. The last perturbative order, g 6 log(l/g) 
has recently been computed using the technique of dimensional reduction [[32|]. 

2.0 




1.5 



1.0 



0.5 



1000 



,1,1,1,1 








Stefan-Boltzmann law 




0(g 6 [ln(l/g)+const.]) 




— 4d lattice data 


J , 1 


i.i.i, 



2.0 



3.0 



4.0 



5.0 



6.0 



Figure 13: The pure gauge EOS in perturbation theory and the comparison with LGT data, from Refs. [ [32[ 
|33l. 



The convergence, for pure gauge theory, is illustrated in Fig. |13J. While the convergence is 
questionable, with fitting the unknown, non-perturbative C(,g 6 term to (pure gauge) LQCD results 



for T > 3T C , good agreement can be found Q33[]. The EOS, at high temperature, is thus reasonably 
well understood. The approach to the Stefan-Boltzmann (SB) limit is rather slow, due to the slow, 
only logarithmic vanishing of g(T) with increasing temperature. 

The high temperature behavior of pure gauge LGT in 2 + 1 dimensions was recently investi- 
gated in detail [34]. In 2 + l-d, the dimensionless ratio g 2 /T serves as the running coupling that 
vanishes at high temperatures. Since the interaction measure (the trace anomaly) vanishes for free 
gluons, one expects 



e-2p g 



T 3 



y /( lQ g — 



(4.1) 



Fig. [14| shows pressure and interaction measure, scaled to the leading perturbative behavior. 
The data nicely confirm that (e — 2p) / (gT) 2 approaches a constant at high temperature. The EOS 
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o 4x32 




Figure 14: The pressure (left) and the interaction measure, scaled to leading perturbative behavior, (right) 
for purge gauge LGT in 2 + 1 dimensions, from Ref. [[m||. 



is, again, well understood at high temperature, increasing our confidence in out understanding of 
the 3 + l-D theory. Due to the faster running of the coupling, the pressure around T ~ 2T C is even 
farther from the SB limit than for 3 + l-d QCD. 



4.3 The EOS around the transition/crossover region 

While the EOS is well understood and modeled at low and high temperatures the non-perturba- 
tive input from lattice calculations is needed to determine the EOS at intermediate temperatures, 
i.e., in the transition/crossover region. It is also needed, as we have seen, to calibrate the high 
temperature description — via the fit of the €?(g 6 ) term — and to check the range of validity of the 
HRG model at low temperatures. 
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Figure 15: Quark mass effects from NLO pertur- 
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of the EOS [33| (see text). 
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Figure 16: Quark mass effects in the lattice EOS 
with staggered quarks. Nf =2 standard staggered 
data are from Refs. [§5j [3(|| and Nf = 2+1 stout- 
link data from Ref. [p] . 



This intermediate temperature range is also the region where the effect of the strange quark, 
with m s /r c ~ ^(1), and even of the charm quark, become visible. This can be seen in Fig. p~5| which 
shows the quark mass effects incorporated as NLO perturbative corrections to a phenomenological 
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modeling of the EOS based on the HRG, pure gauge lattice results and perturbation theory [[33|], and 
Fig. 16 , which compares lattice results for staggered quarks with Nf • = 2 [ 35[ 36 ] and Nf = 2+1 [§] 
flavors 1 . The effect of the strange quark becomes clearly noticeable around r« 1.5!T C and the effect 
of the charm quark above T ~ 3T C . 
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Figure 17: The interaction measure (left) and the energy density (right) from 2+1 flavor simulations with 
asqtad fermions along two lines of constant physics [p7|]. 

There have been two recent computations of the EOS in full QCD, on the path to a determina- 
tion in the continuum limit, both using N t =4 and 6 lattices. The first, using stout-link improved 
fermions and the exact RHMC algorithm, appeared in Ref. [EJ. Since stout-link fermions do not 
have an improved high temperature behavior, the authors of Ref. [^] applied an a-posteriori tree- 
level correction factor c cont /c^ t = 0.571 and 0.663 for N t = 4 and 6 to correct for the tree level 



(a ) lattice effects (see difference in the data presented in Fig. 16 and Fig. 18 (right)). 




Figure 18: The pressure with asqtad fermions, from Ref. [[37]], (left) and with stout-link fermions, corrected 
for tree-level lattice effects, from Ref. [^|], (right) for full QCD simulations. 

The second, by the MILC collaboration, described in more detail in the contribution by L. Lev- 
kova [37], uses asqtad fermions. For the most part, configurations generated earlier with the R- 

1 Incidentally, this figure also indicates that the stout-link improvement has negligible impact for the EOS. 
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algorithm ^ were analyzed. The major effort, over the last six months, concerned an estimate of, 
and correction for, the step-size errors induced by the R-algorithm as described in more detail in 
the contribution by L. Levkova [[3^] . 

The preliminary, step-size corrected results along two different lines of constant physics, with 
mi = 0.2m s and mi = 0.lm s , and two N t 's, 4 and 6 for the latter case, for the interaction measure and 



the energy density are shown in Fig. |17[ Overall, light quark mass and lattice spacing effects appear 
to be fairly small. The pressure along these lines of constant physics is compared with the result 
from stout-link fermions of Ref. [||] in Fig. 18. With appropriate caveats, the two determinations 
agree reasonably well. 

A reliable continuum extrapolation of the EOS for full QCD should become possible in the 
not too distant future. Certainly, simulations with N t = 8 are needed, and a repeat of the asqtad 
fermion simulations for N t =4 and 6 with the exact RHMC algorithm would be desirable. 

5. Other results 



1.2 



0.6 



0.4 




Figure 19: The space-like string tension for pure gauge theory and for 2+1 flavor QCD, from Ref. [|28|]. 

At high temperature, one expects dimensional reduction for QCD to work well, especially 
in the chromomagnetic sector. The dimensionful 3-d gauge coupling is given by g 2 = g 2 (T)T, 
and the fermions, to leading order, affect only the running of g(T) - higher order corrections can 



be included, see Ref. [42]. This proposition was studied recently by measuring the spatial string 
tension in 2 + 1 flavor QCD with p4 fermions and comparing to pure gauge results, as shown in 
Fig. 19. The comparison can be quantified by fitting yfc = csg 2 = cj,g 2 {T)T, with C3 = 0.587(41) 
for 2 + 1 flavor QCD, as compared to C3 = 0.566(13) for 4-d SU(3) pure gauge theory p3| ] and 
C3 = 0.553(1) for the 3-d theory fl44|]. Similarly, for 2-flavor QCD with clover fermions, Ukita 
reported C3 = 0.54(6) [|3|], compatible with these results. 

Several other interesting results in finite temperature lattice gauge theory appeared during 
the last year or were presented at this conference. These include a study of the dynamics of the 



phase transition in pure SU(3) lattice gauge theory [p8|], indications of a non-trivial fixed point in 
finite temperature U(l) lattice gauge theory [t39|], a suggestion that the center vortex field is the 
field that becomes massless at the deconfinement transition of pure gauge SU(2) theory [ffij], and 
impressively accurate results for two-color strong coupling QCD with staggered fermions at finite 



temperature, and finite chemical potential, with a novel cluster algorithm [41] 
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6. Conclusions 



During the last year simulations with improved staggered fermions have produced convincing 
evidence that in nature, at zero baryon density, there is a crossover from the hadronic phase to the 
quark-gluon plasma phase, rather than a genuine phase transition [17]. A new determination found 
a crossover temperature of 192(8) MeV [27], some 10-15% larger than earlier estimates. Both 
results need confirmation from other groups and, in particular, from other fermion discretizations. 

Progress has been made in the computation of the EOS for full QCD at intermediate tem- 
peratures, around the crossover region, where model and perturbative calculations are unreliable. 
A 7 , = 8 results, though, are needed for controlled continuum extrapolations. 

Much remains to be done in understanding the quark-gluon plasma phase in the vicinity of the 
crossover. For example, determinations of transport coefficients and response functions would be 
very desirable. 

Acknowledgments: Thanks to all who provided information and figures for this review, and thanks 
to my colleagues from the MILC collaboration and to Frithjof Karsch for numerous discussions on 
the subject covered here. 
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